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PEAK SWITCH DETECTOR FOR TRANSIT TIME ULTRASONIC METERS 

CROSS-REFERENCE TO RELATED APPLICATIONS 
Not Applicable. 

STATEMENT REGARDING FEDERALLY SPONSORED 
RESEARCH OR DEVELOPMENT 

Not Applicable. 

BACKGROUND OF THE INVENTION 

Field of the Invention 

[001] A disclosed embodiment of the invention relates generally to the detection of errors in 
ultrasonic transit time measurements. Even more particularly, a disclosed embodiment of the 
invention relates to the identification of mistakes in peak selection for the ultrasonic waveforms. 
Description of the Related Art 

[002] After a hydrocarbon such as natural gas has been removed from the ground, the gas 
stream is commonly transported fi:om place to place via pipeUnes. As is appreciated by those of 
skill in tiae art, it is desirable to know witii accuracy the amount of gas in the gas stream. 
Particular accuracy for gas flow measurements is demanded when gas (and any accompanying 
Uquid) is changing hands, or "custody." Even where custody ti^sfer is not taking place, 
however, measurement accuracy is desirable. 

[003] Gas flow meters have been developed to determine how much gas is flowing through the 
pipeline. An orifice meter is one established meter to measure tiie amount of gas flow. Certain 
drawbacks with this meter existed, however. More recently, another type of meter to measure 
gas flow was developed. This more recently developed meter is called an ultirasonic flow meter. 
[004] Figure lA shows an ultirasonic meter suitable for measuring gas flow. Spoolpiece 100, 



suitable for placement between sections of gas pipeline, has a predetermined size and thus 
defines a measurement section. Alternately, a meter may be designed to attach to a pipeline 
section by, for example, hot tapping. As used herein, the term "pipeline" when used in reference 
to an ultrasonic meter may be referring also to the spoolpiece or other appropriate housing across 
which ultrasonic signals are being sent. A pair of transducers 120 and 130, and their respective 
housings 125 and 135, are located along the length of spoolpiece 100. A path 110, sometimes 
referred to as a "chord" exists between transducers 120 and 130 at an angle G to a centerline 105. 
The position of transducers 120 and 130 may be defined by this angle, or may be defined by a 
first length L measured between transducers 120 and 130, a second length X corresponding to the 
axial distance between points 140 and 145, and a third length D corresponding to the pipe 
diameter. Distances D, X and L are precisely determined during meter fabrication. Points 140 
and 145 define the locations where acoustic signals generated by transducers 120 and 130 enter 
and leave gas flowing through the spoolpiece 100 (i.e. the entrance to the spoolpiece bore). In 
most instances, meter transducers such as 120 and 130 are placed a specific distance from points 
140 and 145, respectively, regardless of meter size {i.e. spoolpiece size). A fluid, typically 
natural gas, flows in a direction 150 with a velocity profile 152. Velocity vectors 153-158 
indicate that the gas velocity through spool piece 100 increases as centerline 105 of spoolpiece 
100 is approached. 

[005] Transducers 120 and 130 are ultrasonic transceivers, meaning that they bodi generate and 
receive ultrasonic signals. "Ultrasonic" in this context refers to frequencies above about 20 
kilohertz. Typically, these signals are generated and received by a piezoelectric element in each 
transducer. To generate an ultrasonic signal, the piezoelectric element is stimulated electrically, 
and it responds by vibrating. This vibration of the piezoelectric element generates an ultrasonic 
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signal that travels across the spoolpiece to the corresponding transducer of the transducer pair. 
Similarly, upon being struck by an ultrasonic signal, the receiving piezoelectric element vibrates 
and generates an electrical signal that is detected, digitized, and analyzed by electronics 
associated with the meter. 

[006] Initially, D ("downstream") transducer 120 generates an ultrasonic signal that is then 
received at, and detected by, U ("upstream") transducer 130. Some time later, U transducer 130 
generates a return ultrasonic signal that is subsequently received at and detected by D transducer 
120. Thus, U and D transducers 130 and 120 play "pitch and catch" with ultrasonic signals 115 
along chordal path 110. During operation, this sequence may occur thousands of times per 
minute. 

[007] The transit time of the ultrasonic wave 115 between transducers U 130 and D 120 depends 
in part upon whether the ultrasonic signal 1 15 is traveling upstream or downstream with respect to 
the flowing gas. The transit time for an ultrasonic signal travehng downstream (i.e. in the same 
direction as the flow) is less than its transit time when traveling upstream (Le, against the flow). In 
particular, the transit time ti, of an ultrasonic signal traveling against the fluid flow and the transit 
time t2 of an ultrasonic signal travelling with the fluid flow may be defined: 

t.=-^ (1) 

c-V- 
L 



c + V— 
L 



(2) 



where. 
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c = speed of sound in the fluid flow; 

V = average axial velocity of the fluid flow over the chordal path in the axial direction; 

L = acoustic path length; 

X = axial component of L within the meter bore; 

ti = transmit time of the ultrasonic signal against the fluid flow; and 

t2 = transit time of the ultrasonic signal with the fluid flow. 
[008] The upstream and downstream transit times can be used to calculate the average velocity 
along the signal path by the equation: 

y=]Lhz}± (3) 
2x tjt2 

with the variables being defined as above. 

[009] The upstream and downstream travel times may also be used to calculate the speed of sound 
in the fluid flow according to the equation: 



c = 



(4) 

2 t,t, 



[0010] To a close approximation, equation (4) may be restated as: 

V=-^ (5) 
2x 



where, 

At = tl-t2 (6) 
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So to a close approximation at low velocities, the velocity v is directly proportional to At. 
[0011] Given the cross-section measurements of the meter carrying the gas, the average velocity 
over the area of the meter bore may be used to find the volume of gas flowing through the meter or 
pipeline 100. 

[0012] Li addition, ultrasonic gas flow meters can have one or more paths. Single-path meters 
typically include a pair of transducers that projects ultrasonic waves over a single path across the 
axis (i.e. center) of spoolpiece 100. In addition to the advantages provided by single-path ultrasonic 
meters, ultrasonic meters having more than one path have other advantages. These advantages 
make multi-path ultrasonic meters desirable for custody transfer applications where accuracy and 
reliability are crucial. 

[0013] Referring now to Figure IB, a multi-path ultrasonic meter is shown. Spool piece 100 
includes four chordal paths A, B, C, and D at varying levels through the gas flow. Each chordal 
path A-D corresponds to two transceivers behaving alternately as a transmitter and receiver. Also 
shown is an electronics module 160, which acquires and processes the data from the four chordal 
paths A-D. This arrangement is described in U.S. Patent 4,646,575, the teachings of which are 
hereby incorporated by reference. Hidden from view in Figure IB are the four pairs of transducers 
that correspond to chordal paths A-D. 

[0014] The precise arrangement of the four pairs of transducers may be more easily understood by 
reference to Figure IC. Four pairs of transducer ports are mounted on spool piece 100. Each of 
Ihese pairs of transducer ports corresponds to a single chordal path of Figure IB. A first pair of 
transducer ports 125 and 135 includes transducers 120 and 130 recessed slightly from the spool 
piece 100. The transducers are mounted at a non-perpendicular angle 9 to centerline 105 of spool 
piece 100. Another pair of transducer ports 165 and 175 including associated transducers is 
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mounted so that its chordal path loosely forms an "X" with respect to the chordal path of transducer 
ports 125 and 135. Similarly, transducer ports 185 and 195 are placed parallel to transducer ports 
165 and 175 but at a different "level" {i.e. a different radial position in the pipe or meter 
spoolpiece). Not explicitly shown in Figure IC is a fourth pair of transducers and transducer ports. 
Taking Figures IB and IC together, the pairs of transducers are arranged such that the upper two 
pairs of transducers corresponding to chords A and B form an X and the lower two pairs of 
transducers corresponding to chords C and D also form an X. 
!f [0015] Referring now to Figure IB, the flow velocity of the gas may be determined at each chord 
R: A-D to obtain chordal flow velocities. To obtain an average flow velocity over the entire pipe, the 
% chordal flow velocities are multiplied by a set of predetermined constants. Such constants are well 
H known and were deteraiined theoretically. 

0 [0016] Because the measurement of gas flow velocity and speed of sound depend on measured 
P: transit time, t, it is important to measure transit time accurately. One method and apparatus for 

1 measuring the time of flight of a signal is disclosed in U. S. Patent 5,983,730, issued November 16, 
1999, entitled " Method and Apparatus for Measuring the Time of Flight of A Signal", which is 
hereby incorporated by reference for all purposes. 

[0017] A difficulty that arises in measuring a time of flight exactly is defining when an ultrasonic 
waveform is received. For example, a detected waveform corresponding to a received ultrasonic 
signal may look like that shown in Figure 2. The precise instant this waveform is deemed to have 
arrived is not altogether clear. One method to define the arrival instant is to define it as a zero 
crossing that follows a predefined voltage threshold value for the waveform. However, signal 
degradation due to pressure fluctuations for example may cause the correct zero crossing to be 
misidentified, as shown in Figure 3 (not to scale). Other methods for identifying arrival time may 



54434.01/1787.12700 



6 



also be used, but each is also subject to measurement error by misidentification of the proper arrival 
time. 

[0018] Further, regardless of the method of measuring arrival time, adaptive methods are known 
that modify the measurement criterion or criteria values during measurement of the ultrasonic 
signals. These methods were developed to be more accurate than other methods of measuring 
arrival time but nonetheless still suffer from misidentification of the ultrasonic signal arrival time, 
leading to measurement errors. Further, because of the adaptive nature of the identification method, 

t the misidentification tends to be "locked in", where the detection algorithm tends to lock in the 

h, mistake so that it happens time after time after time. 

S [0019] Although the misidentification of the arrival time of ultrasonic signals has long been 
^: known, previous approaches to identifying peak switching errors are inadequate. For example, 
p ultrasonic meters generally provide a speed of sound reading for the gas flow through the meters. 
R This speed of sound calculation varies with error in the peak selection. However, the variation in 
i the speed of sound measurement is small and may be difficult to identify directly. In addition, the 

variation tends to get smaller as the meter size increases, making uniform diagnostics nearly 

impossible. 

[0020] There therefore exists a need for a method or ultrasonic meter that eliminates 
misidentification of the arrival time of utosonic signals. Ideally, such a method or meter would be 
compatible with existing meters. It would also be ideal if such a metiaod or meter were inexpensive 
to implement. 

SUMMARY OF THE INVENTION 
A method to detect a peak selection error in a waveform includes measuring a first average 
transit time for one or more ultrasonic signals along a first path in a meter, measuring a second 
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average transit time for one or more ultrasonic signals along a second path in the meter (where the 
first and second paths are of different lengths) and determining whether there exists a peak selection 
eiTor based on the length of the first path, the length of the second path, the first average transit 
time, and the second average transit time. 

BRIEF DESCRIPTION OF THE DRAWD^GS 
[0021] For a more detailed description of the prefenred embodiment of the present invention, 
reference will now be made to the accompanying drawings, wherein: 

Figure lA is a cut-away top view of an ultrasonic gas flow meter; 

Figure IB is an end view of a spoolpiece including chordal paths A-D; 

Figure IC is a top view of a spoolpiece housing transducer pairs; 

Figure 2 is a first exemplary received ultrasonic waveform; 

Figure 3 is a second exemplary received ultrasonic waveform; 

Figure 4 is an end view of a four-chord ultrasonic meter showing angular relationships 

between chords; 

Figure 5 is flow chart of one method according to the invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
[0022] A first embodiment of the invention is a method and variants thereof for detecting time of 
arrival errors for ultrasonic signals. A second embodiment is electronics or hardware implementing 
such methods. For example, a processor or microprocessor associated with an ultrasonic meter as 
shown in Figure 1 might run a computer program embodying the disclosed method. Alternately, 
electronics might analyze data provided by an ultrasonic meter after the fact to determine whether 
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an error was made in determining the measured time-of-flight. 

[0023] For a chord A of known length La, it is known that an ultrasonic wave traveling at the speed 
of sound "c" through a homogeneous medium at zero flow in the meter traverses the length of the 
chord La in time tA. tA is also therefore the average transit time of the ultrasonic signal along chord 
A as fluid moves through the meter. tA may not be found, however, by simply averaging the 
upstream and downstream transit times. Instead, .the value of tA may be found algebraically by the 
equation: 

it follows that: 



L 

c = 

t 



(8) 



A 



This is just as true for a second chord B, such that 
L 



c = — 
t 



B- (9) 

B 



[0024] For various reasons, however, the measured gross transit time is not exactly the actual 
transit time of the signal. One reason, for example, that the two times differ is the delay time 
inherent in the electronics associated with each transducer. 
[0025] If total measured time T is defined as: 

T = t + T (10) 

where, 

T = measured or gross transit time; 
t = actual transit time; and 
T = delay time. 
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Then where the delay times are the same for chords A and B, it is known from equation (8) that: 



c = ■ 



(11) 



[0026] It follows from knowledge that the speed of sound for the medium is the same at both chord 
"A" and chord "B" that: 

La(Tb-t) = Lb(Ta-t) (12) 

and 

1^ X= ^bTa-L^T3 (13) 



AL is defined as: 

AL = Lb-La (1^^) 
and it follows that: 



t — 



(15) 



AL AL 
with the variables being defined as above. 

[0027] Of course the transducer delay time for chord A, ta, and the transducer delay time for chord 
B, TB, are not the same. However, these delay times are routinely measured for each pair of 
transducer at the manufacturing stage before the transducers are sent into the field. Since ta and tb 
are known, it is also well known and common practice to cahbrate each meter to factor out 
transducer delay times for each ultrasonic signal. Effectively, xa and tb are then equal to zero and 
therefore the same. The only remaining component to transducer delay time is a misidentification 
of the arrival times for the ultrasonic signals as by peak selection errors. Since the measured transit 
time T is defined as the actual transit time, t, plus delay time, t, actual transit time can be 
substituted for measured transit time T where there is no peak selection error to result in: 
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= 0 



(16) 



[0028] This equation can then be used as a diagnostic to establish whether an error exists in the 
peak selection. It is equation 16 that has general applicabiUty to a broad range of ultrasonic meters 
and signal arrival time identification methods. 
[0029] A variable r] , may then be established: 



La = length of chord A; 
Lb = length of chord B; 

tA = average transit time of ulttasonic signals traveling along chord A (corrected average 
measured transit time); 

tB = average transit time of ultrasonic signals traveling along chord B (corrected average 
measured transit time); and 

AL = Lb-La (where chord B is longer than chord A). 
[0030] If there is a misidentified peak, rj ?t 0. 

[0031] Let us examine the case where the "actual" peak in the received ultrasonic signal has been 
misidentified in chord A. Referring again to Figure 3, an example ultrasonic wave form 300 is 
shown (not to scale). The "true" arrival time 320 is a zero crossing on the waveform as the 
waveform goes from positive to negative polarity. Of course, a processor or the like may 
mistakenly identify either early zero crossing 310 or late zero crossing 330 (or some other zero 
crossing) as the true or actual arrival time of the waveform. If the ultrasonic meter is operating at 
125 kHz, for example, the period of the waveform is 8 jisec (microseconds). Therefore, for this 
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(17) 



where. 



frequency, any peak selection error is on the order of 8 jisec. Obviously, a peak selection error for 
another frequency will have a different value, depending on the period of the waveform. In any 
event, for an error in tA, 

Tl=^(tA+te)-^ (18) 
' AL AL 

where, 

r\ = error indicator 

La, Lb = lengths of chords A, B; 

tA, tB = average transit times across chords A, B at zero flow; 
AL = difference in the lengths of chords A and B; and 
te = error in transit time. 
This is equivalent to: 

' AL AL AL ' 
From equation (16), it is known that 

AL AL 

such that Ti = ^^t, (20) 
AL 

[0032] Based on the value of t], it can be inferred whether there has been misidentification of the 
arrival time. These equations may be adapted to any particular meter. For example, referring to 
Figure 4, there is shown a schematic end-view of an ultrasonic meter. The meter housing 400 
includes chord A, chord B, chord C, and chord D. A central or longitudinal axis is represented at 
the center of the housing by a center point 410. A vertical axis 420 intersects center point 410. A 
line drawn from the center point at an angle of 36° from the vertical axis 420 intersects both chord 
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A and the meter housing 400. A line drawn from the center point 410 at an angle of 72° from the 
vertical axis 420 intersects both chord B and the meter housing 400. Assuming an ideal meter that 
has the transducer with no setback from bore, according to this: 

Lb= V2Dsin(72);and (21) 

La= V2Dsin(36) 

where, 

La= length of chord A; 
Lb= length of chord B; and 
D = meter diameter 
Therefore: 

= 1.618034 (22) 

La 

— may be expressed as La/(Lb-La). 
AL 

Multiplying by results in . Substituting 1.618034 for Lb/La results in: 

1/L^ (Lg /L^ —1) 

= 1.618 (23) 

AL 

Similarly, — may be expressed as Lb/(Lb-La). 
AL 



1/L. . . Lg/L^ 



Multiplying by — ^ results in _ 
1/ L^ C^b'^a ^) 

Substituting 1.618 for Lb/La results in: 



^ = 2.618 (24) 
AL 
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[0033] Obviously, where there is a different angular (or bounce path) relationship between or 
among the chords, the value of these ratios will differ. 

[0034] Assuming no effect from any recess of the transducers from the perfect circle of the 
theoretical meter, the above values for La/AL and Lb/AL are independent of meter size. This 
assumption is largely accurate, with the ratios La/AL and Lb/AL varying with meter size due to 
recess of the transducers back from the spoolpiece-housing. This variation is greater the smaller the 
meter. 

[0035] Because t was defined as zero above as being calibrated or factored out from the meter, and 
because Lb/AL is equal to 2.618 as shown above with respect to equation (20), the equation 
simplifies to: 

Ti = 2.618te (25) 

[0036] If the misidentified peak is one zero crossing away from the "true" zero crossing (the most 
common occurrence), the value for t| is plus or minus 21 ]Lisec, depending on whether misidentified 
zero crossing is earher in time than the true zero crossing (in which case the value for r| is negative) 
or is later in time than the true or correct zero crossing (in which case the value for j] is positive). 
More generally, 

ti = 2.618Ea- 1.618Eb (26) 

where, 

Ea = the error in tA*, and 
Eb = the error in tB. 

If Eb = 0 and Ea = +/-8ns, then r] - +/-21|a5; 
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ff Ea = 0 and Eb = +/-8ns, then t] = +/-13^s; 

If Eb = Ea = +/-8HS, tlien ^ = +/-8ixs. 
[0037] With normal manufacturing tolerances, the errors Eb and Ea are on the order of +/-l|xs, 
which implies an y\ of about +/- 3|is. So long as the eiror from manufacturing tolerances, 
multipUed by the maximum UAL value {e.g. 2.618) is less than one period of the ultrasonic 
waveform, any error measurement less than some predetermined value can be discarded as mere 
manufacturing tolerance. For example, the predetermined value may be any value less than one 
period such as in this case 4^sec. Thus the ^ test can separate peak switching errors from 
manufacturing tolerances. 

[0038] The above derivation is based on an assumption of zero flow through the meter. Of course, 
in a real-worid application, the assumption of a homogeneous medium is not necessarily true at 
very low velocities, due to stratification and convection. At high flow rates, turbulence and 
pressure surges disturb the homogeneity of the medium. Consequendy, use of these exact 
equations may present challenges that would best be avoided. The preferred embodiment 
implements the underiying theory of this method by use of the measured speed of sound, c. 
[0039] It is preferred that tA be based on the average of a batch of upstream and downstream 
measurements for chord A. Similarly, Ib should be based on the average of a batch of upstream and 
downstream transit time measurements for chord B. However, we akeady know from equation (7) 
that tA = La/ca and ts = Lb/cb. It follows that: 



LbL a(Cb ^a) (27) 



ALc^Cb 



where, 

r| = error indicator 
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La, Lb = lengths of chords A and B; 

ca, cb = values for speed of sound measured by chords A and B; and 
AL = difference in the lengths of chords A and B. 

[0040] This calculation presents an additional advantage. Of course, ultimately this computation is 
based on the same variables as the earlier equations. But because a standard ultrasonic meter such 
as that sold by the assignee already calculates speed of sound for each chord, a value for t] may be 
easily computed based on already known or computed information. 

[0041] Figure 5 is a flow chart of a method according to one embodiment of the invention. A 
method to detect a peak selection error in a waveform includes preprocessing of the wave form (e.g. 
filtering and digitizing the waveform) at step 500, measuring a first average transit time for one or 
more ultrasonic signals along a first path in a meter at step 510, measuring a second average transit 
time for one or more ultrasonic signals along a second path in the meter at step 520, where the first 
and second paths are of different lengths, and determining whether there exists a peak selection 
error based on the length of the first path, the length of the second path, the first average transit 
time, and the second average transit time at step 530. Additional steps may also be added. For 
example, at step 540 the measurement error may be corrected for by knowledge of the magnitude 
and direction or polarity of the error either by correction after the fact or by "on the fly" correction. 
If the value of r| is positive, the error occurs in the later in time direction. If the value of r| is 
negative, the error occurs in the earlier in time direction. 

[0042] In a four-chord meter, this r\ comparison test can also be made between chords A and C, D 
and C, and D and B because each of these pairs of chords is comprised of different-length chords. 
The four values of r\ can then be used to pinpoint and correct the peak switching error. 
[0043] Another application is to check that the time-of-flight measurement is correct. For example, 
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the check on time of flight measurement may be used where a zero crossing method of measuring 
the arrival time of an ultrasonic signal is not used. 

[0044] It should be noted that the invention is applicable not only to a four-chord ultrasonic meter 
as shown above but also to other meter designs including bounce-path ultrasonic meters (so long as 
those meters have at least two chords of differing lengths). The use of nomenclature such as La and 
Lb should not be construed to hmit the invention to a particular set of chords. The invention may 
be used with any chords in an ultrasonic meter. It should also be noted that the invention is not 
limited to zero-crossing points, and could also be appUed to peak selection points or any other 
points on the waveform where the picking errors are relatively constant. The invention can also be 
used to verify accurate time measurement within some tolerance regardless of the method of 
measuring time of arrival (although the indication may be weaker or stronger than with zero 
crossing methods). 

[0045] While preferred embodiments of this invention have been shown and described, 
modifications thereof can be made by one skilled in the art without departing firom the spirit or 
teaching of this invention. The embodiments described herein are exemplary only and are not 
Umiting. For example, the invention applies equally well to digitized signals as analog signals. 
Many variations and modifications of the system and apparatus are possible and are within the 
scope of the invention. For example, the invention apphes to any ultrasonic data originating from 
two or more chords having differing lengths and is not Umited to the disclosed four chord meter. 
Accordingly, the scope of protection is not limited to the embodiments described herein, but is only 
limited by the claims which follow, the scope of which shall include all equivalents of the subject 
matter of the claims. 
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